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ABSTRACT 

The  purpose  of  this  research  was  to  show,  both  experiment- 
ally and  analytically,  that  the  stress  concentration  around 
holes  in  laminated  fibrous  composites  is  not  adequately  de- 
scribed by  the  theoretical  solution  for  homogeneous  orthotropic 
plates,  and  that  it  is  a  function  of  hole  size.   It  was  shown 
that  the  gross  laminate  properties,  determined  analytically  from 
individual  lamina  properties,  make  it  impossible  to  express  the 
proper  boundary  conditions  at  a  free  edge.   A  new  expression 
for  the  modulus  of  elasticity  tangential  to  a  free  boundary 
was  developed.   Thus,  the  problem  became  similar  to  that  of  a 
circular  hole  strengthened  by  an  elastic  ring.   Furthermore, 
the  effective  width  of  such  a  ring  was  found  to  depend  upon 
hole  size.   Thus,  the  stress  concentration  factor  was  shown  to 
be  a  function  of  hole  size,  which  was  not  predicted  by  homogene- 
ous orthotropic  theory.   These  results  were  very  satisfactorily 
verified  by  experiment. 
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n  -  material  parameter  given  by  eqn  (II. E. 2) 

X   -  -  transverse  distance  from  center  of  hole  (II.E.l) 
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Y  -  shear  strain 

a  -  lamina  orientation  angle  from  longitudinal  axis 

8  -  angle  from  longitudinal  axis 

P  -  load 

"h  -  thickness 

SUBSCRIPTS 

1  -  longitudinal  axis  of  composite  plate 

2  -  transverse  axis  of  composite  plate 

45  -  45°  between  1  and  2  directions 

l  -  longitudinal  (fiber)  axis  of  a  lamina 
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t  -  transverse  axis  of  a  lamina 

n  -  normal  to  plate  boundary 

s  -  tangential  to  plate  boundary 

6  -  arbitrary  angle  from  longitudinal  axis 

SUPERSCRIPTS 

c  -  composite  plate  (includes  photoelastic  coating) 

c-P  -  composite  plate  less  photoelastic  coating 

P  -  photoelastic  coating 

(k)  -  denotes  K    lamina 


I.   INTRODUCTION 

A.  BACKGROUND 

Because  of  their  high  strength-to-weight  ratios,  laminated 
fibrous  composites  are  expected  to  be  utilized  extensively  for 
aerospace  vehicles  in  the  future.   At  the  present  time,  how- 
ever, the  realization  of  the  maximum  potential  of  composites 
is  severely  restricted  by  the  limited  techniques  available  for 
their  analysis.   An  aerospace  vehicle  is  necessarily  a  very 
complex  structure  which  is  difficult  to  analyze  even  when  com- 
posed of  homogeneous  isotropic  materials.   The  introduction  of 
inhomogeneous  anisotropic  composites  greatly  magnifies  the 
problem.   Lekhnitskii  presents  a  solution  for  stresses  around 
elliptic  holes  in  homogeneous  orthotropic  plates,  but  the  ap- 
plicability to  nonhomogeneous  composites  is  questionable. 

B.  OBJECTIVES 

The  first  objective  of  this  research  was  to  show  experi- 
mentally that  the  homogeneous  orthotropic  solution  for  stress 
concentration  around  a  circular  hole  is  not  accurate  for 
laminated  fibrous  composites.   The  second  objective  was  to 
determine  whether  hole  size  has  an  effect  on  the   stress 
concentration  in  an  infinite  plate.   The  third  objective  was 
to  analyze  the  effects  of  inhomogeneity  on  stress  concentra- 
tion.  Specifically,  the  boundary  conditions  on  a  free  edge 
in  an  inhomogeneous  laminate  were  analyzed  in  general  and 
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in  particular  for  the  composite  tested  experimentally  in  this 
research. 

C.   APPROACH 

As  stated  above,  the  experimental  objectives  required  . 
the  comparison  of  stress  fields  for  different  hole  sizes  with 
the  homogeneous  orthotropic  solution.   A  (0°,-45°)s  glass- 
epoxy  composite  was  used.   Four  different  hole  sizes  were 
tested  by  reflection  photoelasticity .   The  analytical  approach 

was  to  examine  the  relationships  between  constituent,  lamina, 

2 

and  laminate  material  properties  given  by  Calcote  .  The  pur- 
pose of  this  analysis  was  to  determine  if  this  attempt  to  re- 
present laminated  fibrous  composites  as  homogeneous  materials 
was  valid  near  free  boundaries. 
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II.   EXPERIMENTAL  ANALYSIS 

A.   SPECIMEN  DESIGN 

1.  Material  Properties 

The  solution  given  by  Lekhnitskii  [3]  for  the  stress 
concentration  around  a  circular  hole  in  an  orthotropic  plate 
is  dependent  upon  the  material  properties.   Therefore,  the 
first  task  was  to  determine  these  properties  experimentally. 
Due  to  the  thickness  of  the  photoelastic  coating  to  be  used 
for  the  stress  concentration  analysis ,  the  entire  glass-epoxy- 
photoelastic  system  was  analyzed  as  a  whole,  rather  than  try- 
ing to  apply  correction  factors.   Thus,  the  material  properties 
were  also  obtained  by  photoelasticity . 

Determination  of  the  material  properties  of  an  ortho- 
tropic  material  requires  tension  specimens  cut  along  the 
principle  1-t  axes,  and  a  third  specimen  cut  at  45°  to  these 
axes [2],   Therefore,  standard  tensile  specimens  as  described 
in  the  Advanced  Composites  Design  Guide  [1]  were  made  up. 

2 .  Stress  Concentration 

Previous  experimental  work  by  this  author  seemed  to 
indicate  that  stress  concentration  in  composites  is  a  func- 
tion of  hole  size.  One  purpose  of  the  present  research  was 
to  attempt  to  verify  this  hypothesis.  Therefore,  hole  size 
was  the  only  independent  variable  considered  in  designing 
specimens.  All  specimens  were  cut  from  the  same  sheet  of 
material.   The  width  of  each  specimen  was  twelve  times  the 
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hole  radius.   The  dimensions  of  all  specimens  are  given  in 
Table  II. A. 1. 

TABLE  II. A. 1 

1200  SERIES  SPECIMENS 

MATERIAL:   GLASS/EP0XY/PS1-C 
G/E  LAYUP:  (0°,  +  45°  )„ 


THICKNESS:  .15" 


PLATE  # 

LENGTH 

WIDTH 

HOLE 

STRESS  ANGLE 

L 

(in) 

W  (in) 

DIAMETER 
(2a) 

(a) 

1200 

9 

1 

- 

0° 

1210 

9 

1 

- 

45° 

1220 

9 

1 

- 

90° 

1203 

18 

3 

1/2 

0° 

1205 

12 

1-1/2 

1/4 

0° 

1206 

7-1/2 

3/4 

1/8 

0° 

B.   SPECIMEN  PREPARATION 

The  G/E  composite  was  obtained  from  3-M  Company  in  a  two 
by  three  foot  sheet.   The  specimens  were  cut  from  this  sheet 
using  a  diamond  edge  saw.   All  cutting  operations  were  water 
cooled.   The  ends  of  each  specimen  were  built  up  by  attaching 
pieces  of  (0°,  90°)  G/E  as  described  in  the  Composites  Design 
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Guide.   The  most  difficult  operation  was  that  of  drilling 
holes  in  the  G/E  material.   The  difficulty  lay  in  preventing 
fiber  pullout  as  the  drill  bit  came  out  the  back  of  the 
material.   This  problem  was  overcome  by  sandwiching  the 
composite  between  aluminum  when  drilling. 

The  PS-1C  photoelastic  coating  was  cut  one-eighth  inch 
oversize  and  the  holes  were  drilled.   After  cementing  the 
coating  to  both  sides  of  the  G/E  composite,  the  holes  were 
reamed  and  the  edges  were  trimmed  by  machining  at  high  speed 

C.  TESTING  PROCEDURE 

The  photoelastic  data  was  taken  using  the  0  30  series  re- 
flection polariscope  by  Photolastic,  Inc.   The  specimens 
were  mounted  in  a  Reihle  10,000  lb.  testing  machine.   In 
order  to  reduce  the  scatter  due  to  viscoelastic  effects, 
specimens  1200,  1210,  and  1220  were  cycled  to  300  lb.  three 
times  before  testing  began.   The  other  specimens  were  held 
at  their  test  load  for  five  minutes  before  taking  data.   The 
tare  values  were  taken  after  this  initial  stretching. 

D.  MATERIAL  PROPERTIES  RESULTS 

The  original  and  reduced  data  are  given  in  tabular  form 
in  the  Appendix.   Note  that  the  fringe  value,  F,  has  been 
divided  by  the  scale  factor  of  47  of  the  linear  compensator. 
Thus,  the  "fringe  orders"  given  (i.e.  NN,  NO,  TNN,  and  TNO) 
are  actually  compensator  readings.   Using  these  data,  the 
stress-strain  curves  for  each  specimen  and  the  Poisson's 
Ratio  vs.  strain  curves  for  specimens  1200  and  1220  were 
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point  plotted  (Figures  II. D. 1-5).   There  is  considerable  data 
scatter,  particularly  in  Poisson's  Ratio.   Although  a  good  deal 
of  effort  was  directed  towards  eliminating  this  scatter,  visco- 
elastic  effects  and  slight  variations  in  load  during  testing 
made  the  task  extremely  difficult.   It  appears  that  the  visco- 
elastic  properties  of  the  material  affect  the  transverse 
strain  to  an  even  greater  degree  than  the  longitudinal  strain. 
In  view  of  these  difficulties,  due  consideration  was  given  to 
the  relationship, 

E1V21=  E2V12  (II.D.l) 

when  the  curves  were  drawn  in  these  figures.   Using  the  equa- 
tion 

1      4      1      1       1 5> 

+  -f^  (II. D. 2) 


G12    E45     El     E2     El 


the  following  material  properties  were  derived  from  the 
curves: 


El  =  1-24   MSI      v12  "  .5 


c 

2  v21    .31 


E0  =  0.78  MSI       v 


G12  =  .51  MSI 


E.   STRESS  CONCENTRATION  RESULTS 

A  point-by-point  photoelastic  analysis  was  done  for  plates 
1203,  1205,  and  1206.   The  points  chosen  were  along  the  X  and 
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FIGURE  I I.D.I.  -  LONGITUDINAL  MODULUS 
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FIGURE  II. D. 2  -  Modulus  at  45  Degrees 
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FIGURE  II. D. 3  -  TRANSVERSE  MODULUS 
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FIGURE  II. D. 4  -  Longitudinal  Poisson's  Ratio 
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FIGURE  II. D. 5  -  Transverse  Poisson's  Ratio 
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Y  axes  as  shown  in  Figure  II.E.l.   The  data  are  tabulated  in 
the  Appendix  as  explained  in  the  previous  paragraph.   The 
longitudinal  strain  variation  along  the  X-axis  from  the  hole 
to  the  edge  of  the  plate  is  plotted  in  Figures  II. E. 2-5. 

The  theoretical  strain  concentration  at  the  hole  is  given 
by  Lekhnitskii  as 

K  =  p(l+n)  (II.E.l) 

where 


n  =  /  2  (=—  -  v,  0)  +  7t- 


■«,    "12<    G12  (II-E-2' 


These  equations  give  a  theoretical  strain  concentration  factor 
of  3.15  for  this  material.   It  is  evident  from  Figure  II. D. 5 
that  the  strain  concentration  is  a  function  of  hole  size  for 
this  laminated  fibrous  composite,  even  though  the  hole-size- 
to-plate-width  ratio  is  the  same.   Also,  the  measured  value 
seems  to  approach  the  theoretical  value  as  hole  size  increases 
Previous  tests  by  this  author  on  plates  of  constant  width  gave 
the  same  results. 

Figure  II. D. 5  also  shows  that  the  higher  the  strain  con- 
centration a_t  the  hole,  the  faster  it  drops  off  with  the 
normalized  distance  from  the  hole.   In  fact,  the  former  is 
actually  lower  at  some  points!   This  phenomenon  is  necessary 
so  that  the  average  value  of  the  stress  concentration  over 
the  entire  X-axis  (including  the  hole)  will  be  equal  to  one. 
Note  also  the  drop  in  strain  concentration  at  the  outside  edge 
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FIGURE  II.E.l  -  Coordinate  System 
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of  each  plate.   This  effect  is  due  to  a  change  in  gross 
laminate  properties  at  a  free  edge,  which  is  discussed  in  de- 
tail in  Section  III. 
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III.   THEORETICAL  ANALYSIS 

The  solution  given  by  Lekhnitskii  for  the  stresses  around 
elliptic  holes  in  orthotropic  plates  assumes  a  homogeneous 
material.   Because  the  stress  distribution  is  a  function  of 
the  material  properties,  it  is  necessary  to  formulate  gross 
laminate  properties  which  approximate  the  elastic  properties 
of  the  composite  by  those  of  a  homogeneous  system.   The  obvi- 
ous question  here  is,  "Are  the  basic  assumptions  necessary 
for  this  formulation  valid  throughout  the  plate,  particularly 
near  a  free  boundary?" 

A.   BASIC  FORMULATION 

The  basic  formulation  will  first  be  derived  in  detail  be- 
cause it  is  the  basis  for  the  discussion  which  follows.   The 
stress-strain  relationship  (referred  to  natural  axes)  for  an 
orthotropic  lamina  is 

{a}  =  [C]  {e}  (III.A.l) 

where 


{e}  =  fel 


f 


et 


It  J 
and 
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[C]  = 


Cll  C12   ° 


C, „  C^^   0 
12   22 


- 


0    0 


'44 


Referred  to  arbitrary  axes,  this  relationship  becomes  [2] 


=  [C] 


(III. A. 2) 


where 


[C]    =  [T]  1    [C]     [T] 


and  [T]  is  the  matrix  of  direction  cosines 


[T]  = 


cos  a 

.  2 

sin  a 


.  2 

sin  a 

2 

cos  a 


2sina    cosa 


-2sina    cosa 


2        2 

-sina  cosa   sina  cosa   cos  a   -sin  a 


(III. A. 3) 


Equilibrium  requires  that  (assuming  unit  width) 

{P}°  =  Z{P}(k)  =  Z[a] (k)h(k) 

(k)  th 

where  h     is  the  thickness  of  the  k    lamina.   Thus, 


{oY 


Z(a}(k>h(k) 


Eh 


(k) 


(III. A. 4) 


Substituting  equation  III. A. 2  into  this  expression  yields 


"  ih(k> 


(III. A. 5) 


32 


At  this  point,  the  assumption  is  made  that  the  strain  field 
in  each  lamina  is  identical  to  that  of  the  composite  as  a 
whole,  that  is 

{£}   =   (e}C   =  (£}(k)  (III. A. 6) 

equals  a  constant  over  the  summation.   Thus,  equation  III. A. 5 
reduces  to 

[C]C  =  E[C]    k (III. A. 7) 

Eh 

B.   FREE  EDGE  BOUNDARY  CONDITIONS 

The  boundary  condition  for  a  free  edge  in  a  homogeneous 
plate  loaded  in-plane  is 

o=0  (III.B.l) 

n 

By  applying  this  boundary  condition  to  a  composite  represented 
by  the  homogeneous  properties  given  by  III. A. 7,  one  is  actually 
only  requiring  that  the  average  stress  across  the  thickness  be 
zero,  rather  than  the  stress  at  each  point  across  the  thick- 
ness of  the  plate.   Applying  equation  III. A. 4  to  III.B.l, 

Tn    (k),  (k) 

n        ho  h 

ac  =  — 2 =  o  (III.B.2) 

n       EhU) 

(k) 
Thus,  a     in  general  need  not  be  zero,  whereas  the  physics 

of  the  problem  require  that  it  must  be  zero.   Therefore,  the 

material  properties  given  by  equation  III. A. 7  are  not  valid 

near  a  free  boundary  in  a  laminated  composite. 


33 


Equation  III.B.2  results  in  the  free-edge  displacement 
shown  in  Figure  III.B.2.   The  free  edge  is  still  a  plane  nor- 
mal to  the  surface.   Thus,  all  strains  and  displacements  are 

(k) 
constant  through  the  thickness.   The  stress  a     in  a  lamina 

is  derived  as  follows.   The  tangential  stress,  a  ,  is  known 

and  the  normal  stress,  a  ,  is  zero.   Thus,  at  the  free  edge 

e(k)   =  eC  =  ac/EC  (III.B.3) 

s       s    s   s 


and 


(k)      c      c   c 

£       =  £     =  -V     £ 

n       n      sn   s 


But 


or 


(k)     ,    En     :  (k)  ,  c   (k)  c. 

a  =       (t )    (e  +v    e  ) 

n        1-v   v         n   sn   s 
ns  sn 


(k)  =   ( n (k)_vc  )£c  (lll.B.4) 

n       1-v   v         sn   sn   s 
ns  sn 


Thus,  the  error  in  the  boundary  condition  is  directly  propor- 
tional to  the  difference  in  Poisson's  ratios  between  the 
composite  and  the  lamina  in  question.   The  error  is  zero  for 
homogeneous  materials. 

The  free  edge  actually  looks  like  Figure  III.B.3.   The 
displacement  must  be  continuous  between  lamina.   Thus,  the 
displacement  will  vary  through  the  thickness  of  a  given 
lamina. 
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Fig.  III.B.l  -  Coordinate  System 


A_ 


Fig.  III.B.2  -  Displacement  for  aC  = 


n 


0 


Fig.  III.B.3  -  Displacement  for   a(k)-  0 

n 
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C.   MATERIAL  PROPERTIES  AT  A  FREE  BOUNDARY 

It  has  been  shown  that  equation  III.B.2  is  not  the  correct 
boundary  condition  for  a  free  boundary  in  a  laminated  compos- 
ite.  The  correct  condition  is 

a(k)  =  0    for  all  k        (III.C.l) 
n 

In  order  to  derive  the  composite  material  properties,  an  addi- 
tional assumption  is  required: 

eo  =  ejk)  -  e*  (III.C.2) 

s    s      s 

is  constant  through  the  thickness.   Note  that 

ej[k)=  £  (III.C.3) 

n     n 


Hooke ' s  Law  for  this  case  is 


0°   =  E^eo  (III.C.4) 

s     s  s 


and 


a(k)=  E(k)£  (iii.c.5) 


s     s   s 


Recalling  equation  III. A. 4 


or 


v   (k),  (k) 
£o    h 

c      S 


EE(k)e  h<k> 
Ece       s   s 


Ih(k> 
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Thus 


EE(k)h(k) 

Zh(k) 


D.   NUMERICAL  RESULTS 

The  next  step  in  the  analysis  was  to  find  numerical  values 

c  c 

for  E  ,  Young's  modulus  at  a  free  boundary  and  Efl/  Young's 

modulus  in  an  arbitrary  direction  far  from  any  boundaries. 

E   may  be  found  by  use  of  equation  III.C.6,  which  requires  that 

(k) 

E     be  known  for  each  lamina.   A  single  lamina  composite 

could  have  been  made  up  and  tested,  but  the  data  scatter  in- 
herent in  the  properties  of  glass-epoxy  composites  cut  from 
different  sheets  would  have  rendered  the  results  virtually 
useless  for  this  analysis.   Fortunately,  this  procedure  was 
not  necessary.   It  was  possible  to  analytically  reduce  the 
gross  laminate  properties  to  individual  lamina  properties  by 
the  procedure  which  follows.   The  experimentally  determined 
properties  of  the  G/E/PS1-C  composite  are  repeated  here  for 
convenience. 

E^   =1.24  MSI  v^2   =  .5 

E^   =   .78  MSI  v^   =  .31 

G12  =   *51  MSI 
The  stiffness  matrix  may  be  computed  from 
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11 


l-v12v21 


22    l-v12v21 


E2V12 


12 


l-v12v21 


C44  "  G12 


(III.D.l) 


which  give 


[C]  = 


1.47 

.462 
0 


.462 
.923 


0 
0 


.51 


(III.D.2) 
2 


The  properties  of  the  PS1-C  photoelastic  coating  are 

EP  =  .462 
VP  =  .36 
and  the  stiffness  matrix  is 

.190         0 
.528         0 

.169         (III.D.3) 


[C] 


P  _ 


The  thickness  of  the  G/E  and  PS1-C  portions  of  the  composite 
will  also  be  required. 

hC  =  .15 

hC"P  =  .0625 

hP  =  .0875 
The  thickness  of  the  PS1-C  includes  the  cement,  which  has  ap- 
proximately the  same  properties.   Recalling  equation  III. A. 7, 
it  is  possible  to  determine  the  properties  of  the  G/E  system 
without  the  photoelastic  coating. 

hC[C]C  =  hP[C]P+hC~P[C]C_P  (III.D.4) 
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yields 


[C] 


c-p  _ 


2.79 

.843 

0 

.843 

1.48 

0 

0 

0 

.987 


The  inverse  of  equations  III.D.l  are  [2] 

„2 


El  ~  Cll    C 


12 


F   =  f     -  — 

^2     22    C 


22 

2 

12 


11 


(III.D.5) 


v 


12 


12    C 


22 


(III.D.6) 


^11 


G12  "  C44 


which  gives  the  properties  of  the  G/E  system. 


E^   P  =   2.31  MSI 


v12  *   =   .57 


E2C~P  =   1.23  MSI 


v21C"P=   .30 


c-p 


12 


.987  MSI 


Recalling  equation  III. A. 7  and  noting  that  the  thickness  of 
each  lamina  in  the  G/E  system  is  the  same: 


lE]c  =  EJCi" 


n 


(III.D.7) 


where  n  is  the  number  of  lamina. 
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Recalling  that 

[C]  =  [T]"1[C] [T]  (III.D.8) 

it  is  possible  to  solve  for  [C]  because  all  lamina  are  the 
same  except  for  orientation!  Equation  III.D.8  may  also  be 
written  [2] 

4  2     2        4 

C,,=C,,cos  a+2 (C, ~+2C ..) sin  acos  a+C-^sin  a 

2     2         4      4 
C,  =  (C,  ,  +C?:?-4C44)  sin  acos  a+C,?(sin  a+cos  a) 

4  2     2        2 

C22=Cllsin  a+2  ^C]  o+2C44^  sin  acos  a+C?:?cos  a     (III.D.9) 

-  2     2  4      4 

C..,  =  (C,  n+C^n-2C,  0-2C.  ,)sin  acos  a+C.,,(sin  a+cos  a) 
44    11   22    12    44  44 

Noting  that  the  C. .'s  do  not  vary  across  the  summation  in 
equation  III.D.7,  one  may  proceed  as  follows.   The  orientation 
of  the  various  lamina  are 


(1)  -  „(6)  _ 


a    =  a 


0 


a(2)  =  a(5)  =  + 


45' 


„«3)  =  a'4'  =  -45" 


and  equation  III.D.7  reduces  to 


c-p 


1/6 


2     1 


4 
2 

-2 


1    -4 
3     4 


(III.D.8) 
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The    inverse    is 


1.5 

.5 

-1.5 

.5 


1 

0 

1 

-1 


■1.5 
.5 

1.5 
.5 


FcAc-p 


(III.D.9) 


which  gives  the  stiffness  matrix  for  a  single  lamina  as 


[C]  = 


4.78 

.161 

0 

.161 

.852 

0 

0 

0 

.305 


(III. D. 10) 


Applying  equation  III.D.6  gives  the  material  properties  of 
the  lamina. 


£„  =  4.75  MSI 


v£t  =  .19 


et  =   .847  MSI 


vt£  =  .03 


GZt=      *305  MSI 


(k) 


In  order  to  apply  equation  III.C.6,  E     must  be  known  for 
each  lamina.   This  quantity  may  be  computed  from 


or 


eZ  20    ci       .  26    ,  ..  rCZ         0    ,  .  2_6 

— m    =  cos    +  —  Bin    +  1/4  [gj-  2v£t]sm  2 


1rv  =  .211cos29  +  1.18sin26  +  .8sin2(26)    (III. D. 11) 


(k) 
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(k)  c 

E      is  computed  for  each  lamina  and  then  E   is  computed  for 
s  s 

equation  III.C.6. 

For  comparison,  E   may  be  computed  from 


——•  =  cos2  e  +  1.59 -sin2  9  +  .358  sin2  (29) 

or 

— ~   =  .806  cos2  0  +  1.28  sin2  6  +  .289  sin2  (29) 

pc 
9  (III. D. 12) 


c      c 
Eq  and  E   are  plotted  as  a  function  of  0  in  Figure  III.D.l. 


E.   EFFECT  ON  STRESS  CONCENTRATION 

As  seen  in  Figure  III.D.l,  Young's  modulus  in  the  y  direc- 
tion at  the  boundary  is  1.17  MSI  compared  to  1.24  MSI  far  from 
all  boundaries.   The  variation  along  the  X-axis  is  something 
like  Figure  III.E.l.   Thus,  the  material  at  the  boundaries  is 
less  stiff  than  the  material  near  the  center  of  the  plate. 
Now,  suppose  that  the  theoretical  stress  concentration  varia- 
tion for  constant  material  properties  is  given  by  curve  A  in 
Figure  III.E.2.   The  variation  in  material  properties  described 
by  Figure  III.E.l  would  tend  to  reduce  the  stress  concentra- 
tion near  the  boundaries,  because  the  material  is  less  stiff 
in  these  areas  (Curve  B) .   This  effect  is  shown  clearly  at 
the  boundaries  of  the  plates  in  Figures  II. D. 2-5.   Thus,  the 
lower  stiffness  near  a  free  boundary  in  a  laminated  composite 
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requires  that  the  stress  concentration  factor  be  smaller  than 
that  given  by  homogeneous  orthotropic  theory. 

The  effect  described  above  can  also  be  shown  to  vary  with 
the  hole  size.   The  modulus  of  elasticity  of  the  laminated 
composite  is  lowest  at  the  free  boundary  and  increases  to  a 
steady  value  at  a  finite  distance  from  the  boundary  as  shown 
in  Figure  III.E.l.   The  actual  function  will  depend  on  compon- 
ent materials  and  layup  of  composite,  but  will  depend  only  on 
the  X  distance  from  the  free  edge  for  a  particular  composite: 

E,  =  E±     (X)  (III.E.l) 

The  stress  concentration  along  the  X-axis,  however,  is  a 
function  of  the  number  of  hole  radii  away  from  the  boundary: 

K   =  K(X/r)  (III.E.2) 

Thus,  in  terms  of  the  normalized  distance,  X/r,  a  plate  with 
a  small  hole  is  affected  more  by  the  difference  in  stiffness 
than  a  plate  with  a  large  hole.   Therefore,  the  smaller  the 
hole,  the  smaller  the  stress  concentration  factor,  even  in  an 
infinite  plate  or  plates  with  a  constant  hole-size-to-plate- 
width  ratio. 

If  Young's  modulus  in  Figure  III.E.l  is  assumed  to  vary 
like  Curve  B,  rather  than  Curve  A,  one  has  the  problem  of 
stress  concentration  around  holes  which  are  strengthened  by 
elastic  rings.   This  problem  has  been  solved  by  Savin  [4]  who 
describes  similar  results.   The  theoretical  stress  concentra- 
tion for  infinite  copper  plates  with  holes  strengthened  by 
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steel  rings  of  various  widths  is  presented  by  Savin  and  this 
set  of  curves  is  very  similar  to  Figure  II. E. 5,  the  experi- 
mental results  for  plates  with  different  size  holes.   Consider 
a  series  of  such  rings,  one  inside  the  other.   A  large  number 
of  thin  rings  with  appropriate  material  properties  could  be 
used  to  produce  the  variation  in  stiffness  described  by  Curve 
A  of  Figure  III.E.l.   It  is  suggested  that  this  approach  be 
further  investigated  in  order  to  obtain  numerical  results. 
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IV.   CONCLUSION 

The  first  objective  of  this  research  was  to  show  experi- 
mentally that  the  homogeneous  orthotropic  solution  for  stress 
concentration  around  a  circular  hole  is  not  accurate  for 
laminated  fibrous  composites.   Referring  to  Figure  III.E.5, 
the  following  stress  concentration  factors  are  obtained: 

TABLE  IV. 1 

Theory  3.15 

Plate  1203  (r  =  .25")  3.07 

Plate  1205  (r  =  .125")  2.45 

Plate  1206  (r  =  .0625")  2.01 

Obviously,  the  theory  is  not  accurate.   The  second  objective 
was  to  determine  whether  hole  size  has  an  effect  .on  the  stress 
concentration  in  an  infinite  plate.   Referring  to  Table  IV. 1 
again,  it  is  obvious  that  the  stress  concentration  factor  de- 
creases with  hole  size. 

The  third  objective  was  to  analyze  the  effects  of  in- 
homogeneity  on  stress  concentration.   It  was  found  that  to  en- 
force the  boundary  condition  at  a  free  edge  required  the 
modification  of  the  gross  laminate  properties  near  the  edge. 
The  lower  stiffness  obtained  analytically  required  that  the 
stress  concentration  factor  be  lower  than  that  predicted  by 
homogeneous  orthotropic  theory.   Furthermore,  it  was  shown 
that  the  stress  concentration  factor  varies  with  hole  size 

47 


because  the  modified  material  properties  affect  a  larger 
normalized  area  for  smaller  holes. 
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APPENDIX    A:     EXPERIMENTAL    DATA 


A-l       MATERIAL    PROPERTIES     TEST    DATA 


MATERIAL    PROPERTIES    TEST 

PLATE    U    1200 

#    OF    RUNS:    4  AREA       0.150 


F    =       40.243 


9    APR    75 


S  IGMA 


NN 


NO 


DELNN       DELMO       EPSY 


EPSX 


NU 


0 

0 

1.1 

-2.7 

100 

666 

22.9 

20.2 

21.8 

22.9 

535 

-3  72 

.737 

200 

1333 

43.0 

46  .6 

41.9 

49.3 

1239 

-396 

.307 

300 

2000 

62.5 

56.6 

61.4 

69.3 

1712 

-758 

.  443 

400 

2  666 

78.7 

84.4 

77.6 

87.  1 

2135 

-988 

.463 

500 

33:>3 

94.5 

106.1 

93.4 

108.8 

2809 

-949 

.333 

0 

0 

0.4 

-2  .7 

100 

666 

22.6 

19.7 

22.2 

22.4 

458 

-434 

.948 

200 

1333 

41.5 

43.2 

41.1 

45.9 

1116 

-5^7 

.481 

300 

2  000 

61.4 

64.9 

61.0 

67.6 

1626 

-829 

.5  10 

400 

2666 

79.8 

85.6 

79.4 

89.3 

2195 

-1000 

.456 

500 

3333 

96.7 

107.5 

96.3 

113.2 

2777 

-1093 

.395 

0 

0 

-1  .4 

-3.0 

100 

666 

20.9 

21  .0 

22.3 

24.0 

551 

-346 

.628 

200 

1333 

40.2 

42.0 

41.6 

45.0 

1342 

-6  31 

.606 

300 

2  00C 

61.1 

6  5.4 

62.5 

6  8.4 

1613 

-901 

.559 

40  0 

2666 

81.2 

90.4 

82.6 

93.4 

2314 

- 1 0 1  0 

.436 

500 

3  333 

98.1 

106.4 

99.5 

109.4 

2630 

-1404 

.540 

0 

0 

2.0 

1.7 

100 

666 

19.6 

10.2 

17.6 

8.5 

-195 

-903 

<ju  >V  -J*  <** 

•T*  '•£•   'C  T 

200 

133  3 

42.5 

45.0 

40.5 

43.3 

984 

-645 

.655 

300 

2000 

60.3 

66.5 

58.8 

64.8 

1545 

-821 

.531 

400 

2666 

80.  1 

87.6 

78.1 

85.9 

2042 

-1100 

.539 

500 

3333 

99.  1 

109.3 

97.  1 

107.6 

258  7 

-1320 

.510 
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MATERIAL    PROPERTIES    TEST 


6    MAY    75 


PLATE    U     1210 

#    CF    RUNS:    4  AREA       0.150 


40. 248 


S  IGMA 


NM 


NO 


DELNN       DELNO       EPSY 


EPSX 


NU 


0 

0 

-0.3 

-0.8 

100 

666 

17.1 

19.4 

17.4 

20.2 

519 

-181 

.349 

200 

1333 

33.0 

38.3 

33.3 

39.  1 

1020 

-319 

.313 

300 

2000 

49.0 

57.3 

49.3 

58.1 

1523 

-460 

.30? 

400 

2666 

66.0 

77.2 

66.3 

78.0 

2040 

-627 

.307 

500 

3333 

81.7 

94.3 

82.0 

95.  1 

2441 

-859 

.352 

0 

0 

1.4 

-3.5 

100 

666 

18.2 

20.4 

16.8 

23.9 

766 

90 

Jt,  ->.      O-  O 

*i*  ir  T-  v 

200 

1333 

34.  1 

40.5 

32.7 

44.0 

1340 

24 

*r-  f-  t  -r 

300 

2000 

51  .0 

60.4 

49.6 

63.9 

1861 

-134 

.07? 

400 

2666 

66.  0 

81  .0 

64.6 

34.5 

25D1 

-98 

.039 

500 

3333 

81  .6 

98.2 

80.2 

101.7 

2911 

-315 

.108 

0 

0 

-1.3 

-1  .0 

100 

666 

20.8 

23.5 

22.1 

24.5 

589 

-?99 

.508 

200 

1333 

3  3.4 

39.4 

39.7 

40.4 

841 

-755 

.89° 

30  0 

2000 

53.4 

56.5 

54.7 

57.5 

1269 

-931 

.734 

400 

2666 

71.0 

77.0 

72.3 

73.0 

1799 

-1110 

.617 

500 

3333 

83.8 

97.0 

85.1 

98.0 

2491 

-933 

.375 

0 

0 

-1.3 

-1  .0 

100 

666 

24.5 

19.5 

25.8 

20.5 

199 

-839 

•v  ou  -a*  »<- 
^  i-  -r  'i' 

200 

1333 

37.7 

41  .0 

39.0 

42.0 

955 

-603 

.625 

30  0 

2000 

55.3 

58.7 

5  6.6 

59.7 

1326 

-951 

.717 

400 

2666 

71.0 

80.1 

72.3 

81.  1 

1936 

-92  3 

.465 

500 

3333 

87.0 

95.1 

88.3 

95.  1 

2247 

-1306 

.581 
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MATERIAL    PROPERTIES    TEST 
PLATE    #    1220 


6    MAY    75 


#    OF 

RUNS: 

4         AREA       0.15 

•  0               F 

=       40. 

248 

P 

SIGMA 

NN 

NO 

DELNN 

DELNO 

EPSY 

EPSX 

NU 

0 

0 

1.0 

-2.0 

100 

666 

29.8 

28.2 

28.8 

33.2 

664 

-495 

.745 

200 

1  333 

55.9 

65.2 

54.9 

67.2 

1847 

-362 

.196 

300 

2000 

82.7 

9^.6 

81.7 

96.6 

2543 

-744 

.29^ 

400 

2666 

109.0 

127.2 

108.0 

129.2 

3453 

-893 

.259 

0 

0 

2.6 

1  .0 

100 

666 

31.0 

29.8 

28.4 

29.8 

595 

-547 

.919 

200 

1333 

55.4 

66.4 

52.8 

65.4 

182  3 

-301 

.165 

300 

2000 

83.8 

96.8 

81.2 

95.8 

2515 

-752 

.2  99 

400 

2666 

110.0 

130.2 

107.4 

129.2 

3477 

-845 

.243 

0 

0 

3.4 

-1  .2 

100 

666 

32.5 

31.7 

29.1 

32.9 

815 

-356 

.4  37 

200 

1333 

59.0 

64.6 

55.6 

6  5.8 

1734 

-503 

.290 

300 

2000 

86.7 

93.3 

83.3 

99.5 

2654 

-693 

.2  63 

400 

2666 

112.8 

128.6 

109.4 

129.8 

3433 

-969 

.282 

0 

0 

2.6 

-1  .6 

100 

666 

33.5 

28.2 

30.9 

29.3 

555 

-688 

^V  J^-  *V  »«» 

1-  "r  'C-  -■» 

200 

1333 

56.4 

64.5 

53.8 

6S.  1 

1825 

-340 

.186 

300 

2000 

86.6 

96.0 

94.0 

97.6 

251  I 

-369 

.346 

400 

2666 

112.7 

129.0 

110.1 

130.6 

3453 

-973 

.2  63 
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A-2       STRESS    CONCENTRATION    TEST    DATA 


STRESS    CONCENTRATION    DATA 
PLATE    n    1203  X-AXIS 

#    OF    RUNS:     2  LOAD       800.0 


F=       36.596 


2  9    MAY    75 


X/P       TNN 


T\'0 


NN 


NO 


EPSY 


EPSX 


FF 

-2.0 

-6.5 

43.5 

43.8 

1096 

-569 

1.00 

0.52 

1.0 

-8.0 

-6.5 

120.4 

■j.   ~>-  ju  «i/ j.  a. 
n-  t-  'i-  -i-  ^  ^ 

3455 

1*  f  *fr   1*  +C  't- 

3.15 

O-  *&,   J,   -J*    -J*    O 

1.3 

-2.0 

-6.5 

75.0 

75.4 

16  77 

-1139 

1.53 

0.68 

1.5 

-2.0 

-6.5 

55.5 

60.5 

1573 

-530 

1.44 

0.  34 

1.3 

-2.0 

-6.5 

45.6 

49.6 

12  37 

-404 

1.22 

0.30 

2.0 

-2.0 

-6.5 

45.6 

46.4 

1161 

-580 

1  .06 

0.50 

2.5 

-2.0 

-6.5 

45.6 

46.4 

1  161 

-580 

1  .06 

0.50 

3.0 

-2.0 

-6.5 

4  5.5 

46.4 

1  161 

-580 

1  .06 

0.50 

4.0 

-2.0 

-6.5 

45.6 

46.4 

1  161 

-580 

1  .06 

0.5  0 

5.0 

-2.0 

-6.5 

45.6 

46.4 

1161 

-580 

1  .06 

0.50 

6.  0 

-8.0 

-6.5 

30.7 

»•-   -  ■-  »<,  -J{     J,  -JU 

i"  -r  *>-  'r  -r  *»*■ 

1041 

3p,  ,,-.  7^  *,s  -,-.  *,-. 

0.95 

a»  y-  jv  «iu  ju  j; 

FF 

-2.0 

-6.5 

43.5 

43.8 

1096 

-569 

1  .00 

0.52 

1.0 

-5.8 

-6.5 

113.2 

«V  5J5  «**  4c  *JU  Up 

3202 

2.92 

%•*>  .V  «l  J*  JL  «V 

*^  *■(-  *■*  S*  *i*  *t- 

1.3 

-2.0 

-6.5 

75.7 

8  7.^6 

2252 

-627 

2.05 

0.28 

1.5 

-2.0 

-6.5 

56.8 

73.1 

2217 

65 

2.02 

-0.03 

1.8 

-2.0 

-6.5 

47.3 

57.7 

1  720 

-84 

1.57 

0.05 

2.0 

-2.0 

-6.5 

47.3 

50.0 

1  29  7 

-506 

1.  18 

0.39 

2.5 

-2.0 

-6.5 

4  7.3 

50.0 

1297 

-506 

1.18 

0.39 

3.0 

-2.0 

-6.5 

47.3 

50.0 

1297 

-505 

1.18 

0.39 

4.0 

-2.0 

-6.5 

47.3 

50.0 

1297 

-506 

1.  18 

0.  39 

5.0 

-2.0 

-6.5 

47.3 

50.0 

1  297 

-506 

1.18 

0.39 

6.0 

-8.0 

-6.5 

3  5.4 

n*  *r*  *■<-  **-  *v  *r 

1167 

-'-  ul  a,  ^-  v.  ,i, 
*»"*  *i*  -v-  t  nr  *r 

1.05 

jO  J-  *'-  :u  Vj  J, 
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STRESS    CONCENTRATION    DATA 
PLATE    #    1205  X-AXIS 

#    OF    RUNS:    2  LOAD       400.0 


F=       40.248 


28    MAY     75 


X/R       TNN 


TNO 


NN 


NO 


EPSY 


EPSX 


FF 

-0.5 

-5.0 

51  .7 

^53.4 

1424 

-676 

1.00 

0.47      . 

1.0 

-2.5 

-5.0 

117.5 

T-   -V  1*  ;i»  *V  -^ 

3551 

****** 

2.49 

-<-  X.    .».  J-  -».  o, 

T  "i"  *i-  *S»  •¥*  -r 

1.5 

-0.5 

-5.0 

63.2 

82.1 

2694 

130 

1.89 

-0.05 

2.0 

-0.5 

-5.0 

49.9 

51.7 

1394 

-633 

0.98 

0.45 

2.5 

-0.5 

-5.0 

49.9 

51.7 

1394 

-633 

0.98 

0.45 

3.0 

-0.5 

-5.0 

49.9 

51.7 

139^ 

-633 

0.98 

0.45 

4.0 

-0.5 

-5.0 

49.9 

51.7 

13  94 

-633 

0.98 

0.45 

5.0 

-0.5 

-5.0 

49.9 

51.7 

1394 

-633 

0.^8 

0.45 

6.0 

-6.5 

-5.0 

36.1 

JU  -JL.  O*  .A.  %U  *Ar 

v  •»•  r  v  *\-  *r 

1260 

■a,  a*  •*-  -j*  **'  j- 
'i*  »•*  *v  t  -."  -i* 

0.88 

*  *  *  *  *  * 

FF 

-0.5 

-5.0 

51.7 

53.4 

1424 

-676 

1.00 

0.47 

1.0 

-5.2 

-5.0 

111.0 

■j-  ^,  »'-  -■-  j.  *.<* 

*(*   "1*  *i*  '1*  -V  "1* 

3438 

U*  -J,  -4,   *U  •'-  «J, 

*r  *r  -v  *r  v  V 

2.41 

J*  *'.  JU  ■**  -J*  J. 

-r*  *r  *r  '/-  "T-  -r 

1.5 

-5.2 

-5.0 

50.3 

76.0 

2636 

382 

1.85 

-0.14 

2.0 

-5.2 

-5.0 

49.1 

51.0 

1  195 

-990 

0.84 

0.83 

2.5 

-5.2 

-5.0 

49.  1 

51.0 

1  195 

-990 

0.84 

0.8  3 

3.0 

-5.2 

-5.0 

49.1 

51.0 

1  195 

-990 

0.84 

0.33 

4.0 

-5.2 

-5.0 

49.  1 

51.0 

1  195 

-990 

0.84 

0.83 

5.0 

-5.2 

-5.0 

49.1 

51.0 

1195 

-990 

0.84 

0.R3 

6.0 

-2.6 

-5.0 

40.2 

****** 

1266 

»<„  **-  »<,  »■-  A  »•- 
1*  *V  T"  ***  If  T- 

0.89 

J^    JL  *'-  "*'  v *  -'- 

-r-  **r  '.  •  i-  *r  -v 
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STRESS    CONCENTRATION    DATA 
PLATE    #     1206  X-AXIS 

#    OF    RUNS:     2  LOAD       200.0 


F=       36.596 


29    MAY    75 


X/R       TNN 


TNC 


\'N 


NO 


EPSY 


EPSX 


K 


FF 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

FF 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 


1.0 
8.0 


1.0 
1.0 
5.8 
1.0 
8.0 
1.0 
1  .0 
1.0 
1.0 
5.8 


-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 
-5.5 


55. 
109. 

46. 

46. 

46. 

46. 

40. 

55. 
112. 

52. 

52. 

52. 

52. 

41. 


61. 

67. 
57. 
57. 
57. 


*j.  >>,  *t,  j.   -■  .'> 
•V  -i*  *i*  nr  nr*  *** 

61.2 

*Ar  -■-  «'--  -'-•  «Ar  »A» 

*^  "V"  *<*  *!*•  ***  nr1 

64.  I 
53.5 
53.  5 
53.5 

»'  o-  ,'-  ^-  j  »V 


15  90 
3148 
2265 
1710 
1  710 
1710 
1235 
1590 
3229 
1855 
1273 
1273 
12  73 
1275 


-481 

.,;  s£  yji  i^.  1;  ,,; 

523 
-31 
-31 
-31 

a,  j-  a*  ^r  J*  .». 

■v  nn  *r  i*-  *r»  1- 

^-431 

*v  -»'  sp  3p  XC  *,£ 

-109 
-6  91 
-691 
-691 

'i-  V  -»-  "*-  *>%  v 


.00 
.98 
.4? 
.06 
.08 
.08 
.78 
.00 
.03 
.17 
.80 
.80 
.80 
.80 


** 


0.30 


-0.73 

0.02 
0.02 
0.02 

:  $  #  £  if  if. 
0.30 

J.    -  .     -^    -J,    *'- 
•  'r  *i»  v   T*  ** 

0.  06 
0.  54 
0.5  4 
0,  5^ 
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STRESS    CONCENTRATION    DATA 
PLATE    #     1203  Y-AXIS 

#    OF    RUNS:     2  LOAD       300.0 


F=       36.596 


29    MAY    75 


Y/P       TNN 


TNO 


NN 


NO 


EPSX 


EPSY 


FF 

-2.0 

-6.5 

43.5 

43.8 

1096 

-569 

1  .00 

0.52 

1.0 

15.5 

-6.5 

68.7 

•4+  *J»  *V  *&*  ***  '^f 

•v  *r  "  *r  *v»  nr 

-1431 

$###<!# 

-1.31 

if  if  A  *  $  if 

1.3 

-2.0 

-6.5 

20.4 

7.0 

-898 

-78 

-0.82 

-11. 51 

1.5 

-2.0 

-6.5 

14.3 

4.0 

-616 

-20 

-0.56 

-30.80 

1.8 

-2.0 

-6.5 

10.2 

9.3 

-25 

420 

-0.02 

0.0  6 

2.0 

-2.0 

-6.5 

16.7 

14.5 

-215 

46  8 

-0.20 

0.46 

2.5 

-2.0 

-6.5 

26.4 

22.0 

-514 

525 

-0.47 

0.98 

3.0 

-2.0 

-6.5 

30.1 

30.6 

-312 

861 

-0.28 

0.36 

4.0 

-2.0 

-6.5 

38.0 

33.5 

-731 

731 

-0.67 

1.00 

5.0 

-2.0 

-6.5 

40.0 

39.1 

-5  70 

966 

-0.52 

0.5  9 

6.0 

-2.0 

-6.5 

40.0 

39.1 

-570 

966 

-0.52 

0.59 

FF 

-2.0 

-6.5 

43.5 

^43.  8 

1096 

-569 

1  .00 

0.52 

1.0 

14.0 

-6.5 

70.5 

"  ',f  ^  T:  V  '■.' 

-1520 

*r  'f-  *r»  <v  i*  -v^ 

-1.39 

-V  *•-  j,  .i -  J,  „,i„ 

T  T  T  T  T  T 

1.3 

-2.0 

-6.5 

24.3 

13.4 

-832 

129 

-0.76 

6.45 

1.5 

-2.0 

-6.5 

14.4 

2.5 

-706 

-106 

-0.64 

-6.66 

1.8 

-2.0 

-6.5 

14.6 

6.4 

-506 

100 

-0.46 

5.06 

2.0 

-2.0 

-6.5 

22.1 

14.0 

-638 

243 

-0.58 

2.63 

2.5 

-2.0 

-6.5 

26.9 

22.5 

-523 

534 

-0.48 

0.98 

3.0 

-2.0 

-6.5 

33.2 

32.5 

-435 

852 

-0.40 

0.51 

4.0 

-2.0 

-6.5 

38.2 

38.3 

-483 

988 

-0.44 

0.49 

5.0 

-2.0 

-6.5 

42.6 

41  .5 

-629 

1002 

-0.57 

0.63 

6.0 

-2.0 

-6.5 

44.5 

43.5 

-658 

1042 

-0.60 

0.63 
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STRFSS    CONCENTRATION    DATA 
PLATF    #    1205  Y-AXIS 

#    OF    RUNS:     2  LOAD       400.0 


F=       40.248 


28    MAY    75 


Y/R       TNN 


TNO 


NN 


NO 


EPSX 


EPSY 


FF 

-0.5 

-5.0 

51.7 

53.4 

142  4 

-676 

1  .00 

0.47 

1.0 

-5.8 

-5.0 

63.4 

*,~  i,c  v  -,;  -,-  -^ 

-2047 

X.     -»f  .«.  J,  -',  .  <- 

**-  *r  *r  *r  *r  -■" 

-1  .44 

■-'  -  <Jr    J-  -Xr     •.*»  ~'  • 

T*  *»*■  ',*   -1*  *T  'C 

1.5 

-0.5 

-5.0 

2.0 

2.1 

227 

328 

0.16 

-0.69 

2.0 

-0.5 

-5.0 

18.5 

22.4 

124 

889 

0.09 

-0.  14 

2.5 

-0.5 

-5.0 

31.3 

28.5 

-537 

74  2 

-0.38 

0.72 

3.0 

-0.5 

-5.0 

35.2 

34.5 

-489 

947 

-3.34 

0.52 

4.0 

-0.5 

-5.0 

38.0 

41.5 

-291 

1257 

-0.20 

0.23 

5.0 

-0.5 

-5.0 

45.  3 

43.3 

-746 

1034 

-0.52 

0.69 

6.0 

-0.5 

-5.0 

46.7 

46.2 

-708 

1191 

-3.50 

0.59 

FF 

-0.5 

-5.0 

51.7 

53.4 

1424 

^-676 

1  .00 

0.47 

1.0 

-0.8 

-5.0 

65.5 

*v  *f-  *fi  "i*  *?■  n* 

-1962 

I'TYT  'i*  •? 

-1.38 

J-r  JU  -J,    .1,  *J,    »*, 

'i*  -.*  *i*  'r-   *ix  -K- 

1.  5 

-0.5 

-5.0 

4.6 

4.7 

175 

380 

0.12 

-0.46 

2.0 

-0.5 

-5.0 

13.0 

22.4 

165 

909 

0.12 

-0.  18 

2.5 

-0.5 

-5.0 

29.8 

28.5 

-416 

802 

-0.29 

0.52 

3.0 

-0.5 

-5.0 

35.7 

34.5 

-529 

927 

-0.37 

0.57 

4.0 

-0.5 

-5.0 

40.3 

41.5 

-476 

1165 

-0.33 

0.41 

5.0 

-0.5 

-5.0 

44.2 

43.3 

-682 

1116 

-0.48 

0.61 

6.0 

-0.5 

-5.0 

46.4 

46.2 

-684 

1203 

-0.48 

0.57 
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STRESS    CONCENTRATION    DATA 
PLATE    H     1206  Y-AXIS 

U    OF    RUNS:     2  LOAD       200.0 


F=       36.596 


29    MAY    7 


Y/R 

TNN 

TNO 

NN 

NO 

Et>SX 

EPSY 

K 

A 

FF 

-  1.0 

-5.5 

55.6 

61.2 

1593 

-481 

1  .00 

0,30 

1.0 

11.2 

-5.5 

55.6 

Sjs  *,»  -,;;,»-,:  -r- 

-1194 

J,  o-  a,  v«-  J-  »v 
v  -t*  V  -^  -i-  *r 

-0.75 

■V  -/  -i*  -r  5i*  -r 

2.0 

-1.0 

-5.5 

19.8 

44.7 

1233 

1994 

0.78 

-0.62 

3.0 

-1  .0 

-5.5 

33  •  5 

38.0 

-137 

1125 

-0.09 

0.12 

4.0 

-1.0 

-5.5 

43.0 

43.1 

-55? 

1057 

-0.35 

0.52 

5.0 

-1.0 

-5.5 

46.0 

49.0 

-44  3 

1271 

-0.28 

0.35 

6.0 

-1.0 

-5.5 

48.0 

49.0 

-594 

1193 

-0.37 

0.50 

FF 

-1.0 

-5.5 

55.6 

61.2 

15  90 

-481 

1  .00 

0.30 

1.0 

11.2 

-5.5 

55.8 

****** 

-1200 

»«,  «•-  JU  *■-  J*  a* 

■V  *P  "r1  *i*  '<"  "V- 

-0.75 

Ok,   V-  -*  J'   <JL.   _t„ 

'l*  O"  *>*  "r  I1  T- 

2.0 

-1.0 

-5.5 

27.7 

34.  7 

106 

1156 

0.07 

-0.09 

3.0 

-1.0 

-5.5 

36.8 

37.5 

-406 

977 

-0.26 

0.42 

4.0 

-1.0 

-5.5 

45.1 

44.5 

-623 

1057 

-0.40 

0.60 

5.0 

-1.0 

-5.5 

46.6 

44.4 

-744 

997 

-0,47 

0.75 

6.0 

-1.0 

-5.5 

49.0 

47.0 

-777 

1052 

-0.49 

0.74 
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APPENDIX  B:   COMPUTER  PROGRAMS 


B-l   MATERIAL  PROPERTIES  PROGRAM 


C 

C 


100 
110 


120 
130 
140 


30 


150 


155 


160 


170 


EQUIVALENCE  ( T I TLE  ,RTB ( 5 )  5 

REAL*4  TITLE ( 15)  , STRA* 15)  , NN,NO, NU,E PSY * 30)  ,SIGMA*30) 

INTEGER*^  PtEPSX 

INTEGER*^  ITB* 12)/12*0/ 

REAL*4  RTB* 23)/23*0.0/ 

ITB( 1 )=! 

I  TBI 2)=0 

ITB( 3)=6 

ITB(4)=6 

ITB(7)=1 

ITB(  10)  =  2 

ITB( 11 ) =1 

INITIAL  ALPHANUMERIC  DATA 
REA0(5, 100) 
FORMAT* 15A4) 
WRITE(6,  110) 
FORMAT ( 'l',l 
RE  AD  (  5  ,  1  0  0  ) 
WRITE (6,120) 
FOPMATC  ,0«  ,  15A4)  * 
READ<5,130)  K1,A,F 
FORMAT*  lOXt 12  ,3X,F!0 
WPITE(6,140)  K1,A,F 
FORMAT*  '0*  OF  }JNS:  '  ,12 

.3 
DC  30  1=1,10 
READ*5,100)  STRA 
WRITE*6,150)  STRA 

100) 

,12  0) 

100) 

,150) 

i  i 


ALPHANUMERIC 
STRA 

STRA 
5A4) 

TITLE 
TITLE 


5,4X,F10.5) 


AREA  ',F6.3,' 


F  = 


F7 


5, 


READ* 
WRITE*6 
RFAD(5, 
WRITE (6 
FORMAT* 
DO  20  J=l, 


STRA 

STRA 
STRA 
STRA 
,1 5A4) 
Kl 


READ(5,155)  K2 
FORMAT*  54XTI2) 
READ(5,100)  STRA 

TARE    VALUES 
READ(5,160)    P,TNN,TNO 
FORMAT* I 10,2F10.4) 
ISIGMA=P/A 
WRITE*6,170)     P 
FORMATl  '0'  ,  14 

DATA    REDUC 

DO    10    I=1,K2 


,  ISIGMA,TNN,TNO 
,1X, 15,2*  IX, F6.  1)  ) 


10 
180 


READ(5,160)     P,N!M»NO 

ISIGMA=P/A 

SIGMA* I  )  =  FLOAT* I  SIGMA) 

DELNN=NN-TNN 

DELNONO-TNO 

IEPSY=F*<1  .5*OFLN0-DELNN) 

EPSY* I )=FLOAT( IEPSY) 

EPSX=F*( 1 ,5*DELNO-2.*DELNN) 

NU=-EPSX/EPSY* I ) 

WRITE* 6, 180)     P, I SIGMA,NN,NO,DELNN,DELNO,  I E PSY , EPS  X , NU 


FORMAT*  ■      •,T4,IX,I5,4(  IX,  =  6. 
IF( J.EQ.K1 )     Wr  ITE*6,110) 

PLOT    STRESS-STRAIN    CURVE 
IF(J.GT.l)     ITS* 11=2 


I) ,2*2X,IS) , 2X,F4.3) 
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20 


IF(J.EQ.Kl)     ITB(1)=3 

ITB(2)=ITB(2)+1 

IF( ITB( 2 ) . BQ.6)     ITB(2)=1 

CALL    0RAWP(K2»EPSYf SIGMA, IT B,RTB) 

CONTINUE 

STOP 

END 
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B-2   STRESS  CONCENTRATION  PROGRAM 


C 

c 


100 
110 

120 
130 
140 

30 
150 

155 

160 


170 


165 


EQUIVALENCE    ( T I TL E , RTB ( 5 ) ) 

REAL*4    TITLE  < 15)  tSTRAC 15) ,NN,NO, K(30)»X(  30)  t A<?0)  t 

INTEGERS    EPSX,FEPSY 

INTEGER**    ITB(  12  )/12*0/ 

REAL*4    RTB(28)/28*0.0/ 

ITB( 1)  =  1 

ITB(2)=3 

ITB(3)=6 

ITB(4)=6 

ITB(7)=1 

1TB( 10) =2 

ITB( ll)=i 

DO    1     ICCUNT=1,3 

INITIAL    ALPHANUMERIC    DATA 
READ(5,100)     STRA 
FORMAT! 15A4) 
WRITE(6,110)     STRA 
FORMAT* ' 1» T15A4) 
READ(5,100)     TITLE 
WRITE(6,120)     TITLE 
FORMAT( •Ol ,15A4) 
READ! 5,  130)     Kl ,P,F 

FORMAT (  10X,  12,8 X,F 10.5, 4X,F1 0.5) 
WRITE<6,140)     K1,P,F 


KF 


FORMAT*  »0#    OF    RUNS:*  ,12,' 

3) 
DC    30    1=1 


LOAD     «, F6.lt1 


F  = 


F7. 


READ(5i 
WRITE (6 

READ(5, 
W  R I T  E  (  6 
READ(5, 
WRITE(6 
FORMAT ( 
DO  20  J  =  l  , 


10 


100) 
t  150) 
100) 
,120) 

ICO) 

,  150) 
i  i 


STRA 

STRA 
STRA 

STRA 
STRA 
STRA 
,15A4) 
Kl 


READ( 5 
FORMAT 
READ(5 
FA 
READI5 
FORMAT 
DELNN= 
DELNO= 
FEPSY= 
EPSX=F 
KF=1.0 
AF=-EP 
WRITE ( 
FORMAT 
F6 
DA 
DO  10 
READ{5 
FORMAT 
DELNN= 
DELNO= 
IEPSY= 
EPSX=F 
K  (  I  )  =-- 1 
A(  I)=- 
IF (NO. 
IEPSY= 
N0=10* 
EPSX=1 
K( I)=I 


,155)    K2 

( 54X, 12) 

,100)     STRA 

R     FIELD    VALUES 

, 160)     TNN, TNOtNNtNO 

( 10X,4F10.4) 

NN-TNN 

NO-TNO 

F*( 1.5*DELN0-DELNN) 

*<  1  ,5*DELNC-2.*DELNN) 


LOAT(FEPSY) 

0)     TNN,TNC,NN,NO,FEPSY,EPSX,<F ,  AF 

F       ',2(F5.i,lX),2(F6.1,lX),2(I6,lX),F5.2,lX, 


XI  I)  , TNN,TNO,NN,NO 
4) 


SX/F 
6,  17 
(  'OF 
.2) 
TA    R 

I  =  1 1 

,  165 
(5F1 
NN-T 
NO-T 
F*(l 

*  (  1  . 5* 
EPSY/F 
EPSX/F 
NE.0.0 
F*DELN 

*  8 

0**8 
EPSYZFLOAT(FEPSY) 


EDUCTION 

K2 

) 

0. 

NN 

NO 

.5 


*DELNO-DELNN) 
DELN0-2.*DELNN) 

LOAT (FEPSY) 
LOAT( IEPSY) 
)     GO    TO    13 
N  /  1  .  3  6 
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180)  X( 
• ,F3.1 
F6.2  ) 
K 1  )  W  R I T  E  (  6 


I)»T 

,  IX, 


NN,TN0, 
2  (  F  b  .  1  , 


110) 

CJPVE 


MI)  =  10* 

10  WRITE (6, 

180  FORMAT ( • 

tlX, 

IFU.EQ. 

PLOT  STRESS-STRAIN 
IF(J.GT.l)  IT3(1)=2 
IF( IC0UNT.NE.3)  GO  TO  2 
IF(J.EQ.Kl)  ITB(1)=3 
2  CONTINUE 

CALL  DRAWP(K2,X,Kf ITBfRTB) 
20  CONTINUE 

ITB(2)=ITB(2)+1 
IF(ITB( 2)  .FD.6) 
1  CONTINUE 
STOP 
END 


NN, 
IX) 


NO, IEPSr,EPSX,K( I ), A( I ) 
,2(F6.1 ,1X) ,2( 16, IX) ,F5.2 


IT8( 2) =1 
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